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The  Escaig  model  for  thermally  activated  cross  slip  in  face  centered  cubic 
(fee)  materials  assumes  that  cross  slip  preferentially  occurs  at  obstacles 
that  produce  large  stress  gradients  on  the  Shockley  partials  of  the  screw 
dislocations.  However,  it  is  unclear  as  to  the  source,  identity  and 
concentration  of  such  obstacles  in  single  phase  fee  materials.  Embedded 
atom  potential,  molecular  statics  simulations  of  screw  character  dislocation 
intersections  with  120°  forest  dislocations  in  fee  Ni  are  described  that 
illustrate  a  mechanism  for  cross  slip  nucleation.  The  simulations  show  how 
such  intersections  readily  produce  cross  slip  nuclei  and  thus  may  be 
preferential  sites  for  cross  slip.  The  energies  of  the  dislocation  intersection 
cores  are  estimated  and  it  is  shown  that  a  partially  cross  slipped  config 
uration  for  the  intersection  is  the  most  stable.  In  addition,  simple  three 
dimensional  dislocation  dynamics  simulations  accounting  for  Shockley 
partials  are  shown  to  qualitatively  reproduce  the  atomistically  determined 
core  structures  for  the  same  dislocation  intersections. 

Keywords:  cross  slip;  dislocation  intersections;  atomistic  simulations; 
embedded  atom  potentials;  dislocation  dynamics  simulations 


1.  Introduction 

There  is  an  increasing  recognition  of  the  need  to  incorporate  physics-based  models 
of  deformation  in  design  of  structural  components.  Whereas  models  for  predicting 
yield  strength  and  creep  behavior  are  beginning  to  incorporate  significant  physics, 
models  for  fatigue  and  ultimate  strength  remain  mostly  empirical.  This  is  because 
strain-hardening  and  fatigue  resistance  are  highly  influenced  by  dislocation 
micromechanisms,  e.g.  cross-slip,  and  including  physics-based  cross-slip  processes 
in  mesoscale  simulations  (e.g.  dislocation  dynamics)  has  been  difficult.  The  early 
work  of  Escaig  remains  the  most  widely  cited  and  used  model  for  cross-slip 
(see  Eigure  1)  [1^];  however,  this  model  poses  several  difficulties  with  respect  to 
quantitative  simulations.  The  model  is  highly  sensitive  to  choice  of  parameters  that 
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(a) 


Figure  1.  Illustration  of  the  aetivation  energy  for  cross  slip  in  Escaig’s  model:  separation 
of  a  Stroh  constriction  into  a  negative  and  a  positive  constriction. 


have  thus  far  been  difficult  to  obtain.  For  example,  the  constriction  width  required 
for  cross-slip  is  unknown  but  significantly  influences  the  energetics  of  the  cross-slip 
process  [5].  This  difficulty  has  led  to  an  ad  hoc  postulate  that  obstacles  always 
exist  in  materials  and  that  they  enable  sufficient  dislocation  core  constriction  under 
stresses,  thereby  ensuring  cross-slip  [1].  Aside  from  being  unsatisfactory,  this  forces 
cross-slip  models  (particularly  for  single  crystals)  to  make  arbitrary  assumptions 
about  dislocation  obstacle  spacings.  Advances  in  atomistic  simulations  make  it 
possible  to  gain  insights  into  the  cross-slip  process  and  may  serve  to  inform 
mesoscale  simulations  to  accurately  capture  the  atomic-level  physics  of  that 
dislocation  process. 

Atomistic  simulation  works  to  date  have  been  limited  to  calculating  constriction 
energies  and  energetics  of  cross-slip  using  a  single  dislocation  in  a  periodic  unit 
cell  [6,7].  These  simulations  helped  to  unravel  some  interesting  and  potentially  useful 
aspects  of  cross-slip.  The  most  important  finding  was  that  there  exists  configurations 
where  one  of  the  pair  of  constrictions  has  a  negative  energy  (negative  constriction) 
and  the  other  has  a  large  positive  energy  (positive  constriction)  [6,7].  Thus, 
spontaneous  nucleation  of  cross-slip  is  possible  so  long  as  only  one  constriction 
is  needed,  as  is  the  case  for  surface  cross-slip  nucleation  [8].  These  previous 
simulations  found  that  the  total  energy  of  a  pair  of  constrictions,  required  for  cross¬ 
slip  within  the  bulk  of  a  solid,  from  a  pre-existing  constriction  (Stroh  constriction, 
see  Figure  1)  is  always  positive.  The  energy  of  a  Stroh  constriction  can  be  shown 
to  be  a  mean  of  the  energies  of  the  negative  and  positive  constrictions  [9]. 
The  simulations  predict  relatively  high  energies  for  cross-slip  ( >  1  eV)  that  thermal 
activation  would  require  substantial  stress  to  activate  at  temperatures  where  cross¬ 
slip  is  observed  experimentally.  This  does  not  address  the  problem  of  having  to 
assume  ad  hoc  obstacles  that  permit  stress  build-up  in  mesoscale  simulations. 

In  the  early  1960s,  Washburn  [10]  suggested  that  cross-slip  in  fee  materials 
preferentially  occurs  at  dislocation  intersections  with  other  forest  dislocations 
residing  on  the  cross-slip  plane  that  form  collinear  locks  and  glide  locks  [10,11]. 
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Table  1.  Reactions  and  locks  that  occur  between  a  1/2[110]  screw  dislocation  and  forest 
dislocations  on  the  (111),  (111)  and  (111)  planes.  The  screw  dislocation  is  considered  to  reside 
completely,  either  on  the  glide  (111)  or  the  cross  slip  (111)  plane. 
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However,  there  has  been  little  to  no  investigation  of  those  ideas  since  that  time. 
In  this  work,  we  expand  on  the  ideas  of  Washburn.  We  consider  screw  dislocation 
intersections  with  forest  dislocations  residing  on  the  two  other  {111}  type  planes, 
other  than  forest  dislocations  residing  on  the  cross-slip  {111}  type  plane.  If  we  take 
the  screw  dislocation  as  having  a  1/2[110]  Burgers  vector  residing  on  the  (111) 
plane  with  the  (111)  plane  being  the  cross-slip  plane,  the  intersections  that  are 
of  importance  for  cross-slip  are  intersections  with  the  (111)  and  (111)  planes.  Table  1 
gives  the  locks  that  are  formed,  as  well  as  the  reaction  product,  for  intersections 
with  different  Burgers  vectors  on  these  two  planes  as  well  as  the  cross-slip  plane. 
In  Table  1,  the  locks  that  are  formed  are  written  down  for  two  cases:  (i)  the  screw 
dislocation  is  completely  residing  on  the  glide  (111)  plane  and  (ii)  the  screw 
dislocation  is  completely  residing  on  the  cross-slip  (111)  plane.  From  Table  1,  it  is 
clear  that,  by  symmetry,  the  intersections  with  the  (111)  plane  are  identical  to  the 
intersections  with  the  (111)  plane,  and  therefore,  one  needs  to  consider  only  the  screw 
dislocation  intersections  with  the  (111)  plane.  Also,  considering  that  the  screw 
dislocation  can  completely  reside  on  the  glide  (111)  or  the  cross-slip  (111)  plane, 
the  intersection  with  the  '/2[101]  Burgers  vector  on  the  (111)  plane  can  be  shown 
by  symmetry  to  be  equivalent  to  the  intersection  with  the  !/2[011]  Burgers  vector 
on  the  (111)  plane.  As  a  result,  one  needs  to  consider  only  the  intersections  with  the 
'/2[l01]  and  the  '/2[IlO]  Burgers  vectors  on  the  (ill)  plane.  For  the  !/2[l01]  Burgers 
vector,  the  intersection  results  in  a  Lomer-Cottrell  barrier  if  the  screw  dislocation 
is  on  the  cross-slip  (111)  plane  or,  a  glide  barrier  if  the  screw  dislocation  is  residing 
on  the  glide  (111)  plane.  For  the  '/2[110]  Burgers  vector,  the  intersection  results  in 
a  Hirth  lock  irrespective  of  whether  the  screw  dislocation  resides  on  the  glide  plane 
or  the  cross-slip  plane.  The  common  direction  (junction  direction)  between  the  screw 
direction  glide  plane,  (111)  and  the  intersecting  (111)  plane  is  [101]  and  the  common 
direction  between  the  screw  direction  cross-slip  plane,  (111)  and  the  (111)  plane 
is  [Oil].  Due  to  the  fact  that  the  junction  direction  as  well  as  the  type  of  lock  that  is 
formed  is  dependent  on  whether  the  screw  dislocation  resides  on  the  glide  plane 
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or  the  cross-slip  plane,  the  energies  of  these  intersections  may  be  expected  to  depend 
strongly  on  the  plane  in  which  the  screw  dislocation  resides.  This  may  enhance  the 
possibility  of  local  cross-slip  at  the  intersection  region. 

In  this  work,  the  possibility  of  cross-slip  at  dislocation-dislocation  intersections 
was  selectively  examined  using  large-scale  atomistic  simulations  that  contain  more 
than  one  dislocation.  A  glide  dislocation  that  intersects  a  pair  of  forest  dislocations 
is  modeled  and  it  is  shown  that  several  conditions  exist  where  cross-slip  can  nucleate. 
This,  we  believe,  has  the  potential  of  allowing  simulations  of  cross-slip  in  larger  scale 
discrete  dislocation  dynamics  simulations  and  provides  a  better  physical  basis  for 
cross-slip  and  perhaps  to  provide  a  realistic  statistical  representation  of  cross-slip 
during  monotonic  or  cyclic  deformation. 

We  have  performed  molecular  statics  simulations  of  screw-character  dislocation 
intersections  ('/2[110]  Burgers  vector)  with  both  the  '/2[101]  and  the  1/2[110]  Burgers 
vectors  on  the  (111)  plane.  We  have  considered  several  line  orientations  for  the 
intersecting  dislocation  over  a  180°  range  (all  orientations  where  the  intersection 
is  expected  to  be  attractive).  For  the  sake  of  clarity,  we  restrict  ourselves  to  results 
only  for  the  screw  (and  near-screw)  dislocation  intersection,  with  the  V2[101]  Burgers 
vector  residing  on  the  (111)  plane  with  the  120°  (Oil)  line  orientation.  Section  2 
describes  the  simulation  technique  and  the  interatomic  potential  used  for  the 
simulations,  and  Section  3  presents  the  core  structures  and  the  energies  obtained 
from  atomistic  simulations  and  dislocation  dynamics  (DD)  simulations  for  these 
specific  intersections.  Finally,  Sections  4  and  5  give  a  discussion  and  summary, 
respectively,  of  the  results. 


2.  Simulation  technique 

The  atomistic  simulations  described  here  employed  the  three-dimensional  (3D) 
parallel  molecular  dynamics  code,  LAAMPS  [12],  developed  at  Sandia  National 
Laboratory.  A  schematic  of  the  simulation  cell  used  in  the  atomistic  simulations  is 
shown  in  Figure  2.  The  simulation  cell  is  a  rectangular  parallelepiped  with  the  x-axis 
oriented  along  [110],  p-axis  along  [112]  and  the  z-axis  along  [111].  The  dimensions 
of  the  simulation  cell  are  62.0  nm  along  the  x-axis  and  31.5  nm  along  both  the  y-  and 
z-axes  corresponding  to  a  simulation  cell  of  5,405,160  atoms.  A  '/2[110]  screw- 
character  dislocation  is  inserted  in  the  middle  of  the  simulation  cell  using  its 
anisotropic  elastic  displacement  field  [13].  Similarly,  two  120°  dislocations,  having 
a  Burgers  vector  of  '/2[101]  and  line  directions  (011)  spaced  30.0 nm  apart 
(at  x=— 15.0  and  15.0nm,  with  y,  z~0),  were  also  introduced  into  the  simulation 
cell  using  their  anisotropic  elastic  displacement  fields.  The  origin  for  the  initial 
anisotropic  displacement  field  of  the  screw  and  120°  dislocations  were  varied 
(ten  different  origins  were  used)  to  obtain  several  different  core  structures  for  the 
screw  dislocation  intersection.  For  simplicity,  fixed  boundary  conditions  were 
applied  along  all  three  directions  and  energy  minimization  was  performed  using  the 
conjugate  gradient  technique.  A  couple  of  simulations  were  performed  using 
periodic  boundary  conditions  along  the  screw  or  ‘x’  direction.  In  these  cases,  the 
two  intersecting  dislocations  were  opposite  in  character  (-1-120°  and  —60°),  such  that 
periodic  boundary  conditions  could  be  applied. 
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Figure  2.  A  schematic  of  the  simulation  cell  used  in  the  atomistic  simulations  described  in  this 
manuscript.  The  axes  dimensions  are  in  units  of  A. 


Table  2.  Lattice  parameter,  Cq,  cohesive  energy,  E^.,  elastic 
constants  Cn,  Cij  and  C44,  and  stacking  fault  energy,  y,  given 
by  the  Ni  Angelo  Moody  Baskes  EAM  potential. 

Property 

Value 

flO 

0.352  nm 

c„ 

2.464  X  10"Nm^2 

C\2 

1.473  X  10"  Nm"^ 

C44 

1.248  X  10"Nm-2 

E, 

4.45  eV 

Y 

0.089  J  m^2 

2.1.  Interatomic  potential 

The  embedded  atom  potential  used  in  the  simulations  was  developed  for  fee  Ni 
by  Angelo  et  al.  [14]  based  on  the  Voter  and  Chen  format.  Table  2  gives  the  lattice 
parameter,  cohesive  energy,  elastic  constants  and  stacking  fault  energy  given  by  the 
potential  for  fee  Ni.  The  stacking  fault  energy,  ~90mJm“^,  lies  inbetween  a  value 
of  60  and  120mJm“^,  the  values  given  by  the  two  Ni  potentials  that  Rao  et  al.  [6] 
used  in  previous  atomistic  simulations  of  bulk  cross-slip. 

2.2.  Depiction  of  core  structures 

In  order  to  illustrate  the  relaxed  intersection  geometries,  we  take  advantage  of  the 
increase  in  atomic  energy  produced  by  the  strain  field  of  the  partial  dislocations. 
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By  plotting  the  atoms  with  assigned  energies  of  greater  than  — 4.42eV  (the 
approximate  energy  of  atoms  in  the  stacking  fault  region)  the  partial  dislocations 
can  be  visualized  easily  even  in  these  large  simulation  cells.  In  order  to  illustrate  the 
cross-slipped-segment  products  of  the  intersections,  the  positions  are  shown  in  a 
[112]  projection  (i.e.  the  x-z  plane)  as  well  as  the  [111]  projection  (i.e  the  x-y  plane). 
In  the  [112]  projection  segments  spread  on  the  initial  (111)  plane  appear  as  a  single 
line  and  cross-slipped  segments  (i.e.  on  a  (111)  plane)  appear  as  a  pair  of  partials 
separated  by  a  stacking  fault.  The  projection  along  the  j-axis  provides  a  quick 
indicator  of  whether  the  screw  dislocation  is  residing  on  the  (111)  glide  plane  or  the 
(111)  cross-slip  plane  or  a  combination  of  both. 


2.3.  Dislocation  dynamics  simulations 

The  core  structures  obtained  for  the  intersection  using  atomistic  simulations  were 
reproduced  using  a  simple  dislocation  dynamics  (DD)  method,  DDLab,  originally 
developed  at  Lawrence  Livermore  National  Laboratory  [15].  DDLab  runs  within 
MATLAB  and  is  intended  for  simulations  having  a  limited  number  of  dislocation 
lines  while  taking  into  account  the  Shockley  partials  [16].  A  'A  [110]]  screw 
dislocation  and  a  120°  intersecting  dislocation  (line  direction  (Oil),  Burgers  vector 
V2[101])  split  into  their  respective  Shockley  partials  (2.0 nm  apart),  70 nm  in  length 
and  fixed  at  both  ends,  were  initially  introduced  into  an  infinite  simulation  cell  and 
allowed  to  relax  under  no  applied  stress.  The  screw  character  dislocation  was  initially 
taken  to  be  either  completely  dissociated  on  the  glide  (111)  plane,  completely 
dissociated  on  the  cross-slip  (1 1 1)  plane  or,  partly  dissociated  on  the  glide  (111)  plane 
and  partly  on  the  cross-slip  (111)  plane.  The  120°  dislocation  was  initially  taken  to  go 
through  the  origin  of  the  simulation  box,  whereas  the  screw  dislocation  was  initially 
taken  to  be  ±3.0  nm  away  from  the  origin  along  the  ‘j’-axis  ([112]).  The  relaxed 
configurations  were  compared  with  results  from  atomistic  simulations. 


3.  Atomistic  core  structures  for  the  intersection 

Figure  3  depicts  [111]  and  [112]  projections  (i.e.  the  ‘x-y’  and  ‘x-z’  planes)  of  one 
of  the  core  structures  obtained  for  the  screw  dislocation  —120°  dislocation 
intersection,  glide  lock  (Gll)  configuration.  In  this  case,  the  screw  dislocation  was 
initially  positioned  a  [112]  unit  along  the  positive  ‘y’-axis  (±0.8  nm),  away  from  the 
origin.  The  original  placement  of  the  threading  and  cutting  dislocations  allows 
the  dislocation  cores  to  interact  and  combine  as  dictated  by  the  geometry  and 
Peach-Koehler  forces  as  generated  thorough  the  interatomic  potentials.  The  screw 
dislocation  fully  resides  on  the  glide  (111)  plane  and  is  constricted  at  both  the 
intersections  (Stroh  constrictions  [17])  around  one  of  the  Shockley  partials  of  the 
120°  '/2[101]  dislocation  residing  on  the  (111)  plane.  Due  to  the  fixed  boundary 
conditions  along  ‘x’  (which  mimic  a  jog  on  the  screw  dislocation  at  the  fixed 
boundaries),  the  screw  dislocation  is  constricted  in  its  glide  plane  at  both  ends  near 
the  boundary,  which  is  equivalent  to  a  full  Stroh  constriction  [17].  Therefore, 
in  order  to  obtain  the  estimates  of  the  energy  for  the  intersection  core  structure 
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Figure  3.  [Ill]  (jc  >>  plane)  and  [112]  (x  z  plane)  projections  of  the  core  structure  (Gil)  for 
a  screw  120°  intersection  in  fee  Ni  from  atomistic  simulations.  Atoms  with  energy^  greater 
than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in  units  of  A. 


Gil  without  any  jogs,  the  energy  of  a  ‘Stroh’  constriction  needs  to  be  subtracted 
from  the  energy  of  this  core  structure. 

Figure  4  depicts  two  of  the  core  structures  obtained  for  the  screw  dislocation 
—  120°  dislocation  intersection,  Lomer-Cottrell  lock  configurations  (LCl  and  LC2), 
that  were  obtained  using  a  slightly  different  elastic  center  for  the  initial  anisotropic 
displacement  field  of  the  screw  and  120°  dislocations,  as  compared  to  Figure  3. 
The  [111]  and  [112]  projections  are  shown  for  each  case.  In  Figure  4a,  the  screw 
dislocation  was  initially  positioned  a  [112]  unit  along  the  negative  ‘/-axis  (4-0.8  nm), 
away  from  the  origin.  In  Figure  4b,  the  screw  dislocation  was  initially  positioned 
a  [112]  unit  along  the  positive  ‘/-axis  (— 0.8nm),  away  from  the  origin.  After 
relaxation.  Figures  4a  and  b,  the  screw  dislocation  fully  resides  on  the  cross-slip 
(111)  plane.  In  Figure  4b,  as  in  Figure  3,  the  screw  dislocation  is  constricted  at  both 
the  intersections  around  one  of  the  Shockley  partials  of  the  120°  '/2[101]  dislocation 
residing  on  the  (111)  plane.  In  Figure  4a,  the  £r/6[112]  Shockley  partial  of  the  screw 
dislocation  combines  with  the  a/6[121]  Shockley  partial  of  the  120°  dislocation 
to  form  an  extended  node  with  a  resulting  short  Burgers  vector  of  the  type  1/6[01 1], 
or  Lomer-Cottrell  barrier.  At  the  other  end  of  the  Lomer-Cottrell  barrier,  the 
screw  dislocation  is  constricted  to  form  one  half  of  a  Stroh  constriction.  Since  the 
interaction  between  the  1/6[011]  Burgers  vector  and  the  other  Shockley  partial 
of  the  120°  dislocation  (a/6[211])  is  weak,  the  120°  dislocation  needs  to  constrict 
significantly  to  attain  a  lower  separation  distance  at  the  extended  node  reaction 
region.  Also,  as  in  Figure  3,  due  to  the  fixed  boundary  conditions  along  ‘x’,  the  screw 
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Figure  4.  [1 1 1]  (x  7  plane)  and  [112]  (x  z  plane)  projections  of  the  core  structures  (LCl  and 
LC2)  for  a  screw  120°  intersection  in  fee  Ni  from  atomistic  simulations.  Atoms  with 
energy  greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in 
units  of  A. 
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Figure  5.  [Ill]  (x  y  plane)  and  [112]  (x  z  plane)  projections  of  the  core  structures  (PCSl 
and  PCS2)  for  a  screw  120°  intersection  in  fee  Ni  from  atomistic  simulations.  Atoms  with 
energy  greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are 
in  units  of  A. 


dislocation  forms  the  equivalent  of  a  ‘Stroh’  constriction  at  its  ends  and  this  energy 
has  to  be  subtracted  to  obtain  the  true  energy  of  intersection  core  structures  LC 1  and 
LC2  without  any  jogs. 

Figure  5  depicts  two  of  the  relaxed  core  structures  for  the  screw  dislocation 
—  120°  dislocation  intersection,  partially  cross-shpped  locks  (PCSl  and  PCS2),  that 
were  obtained  using  a  slightly  different  elastic  center  for  the  initial  anisotropic 


8 

Distribution  Statement  A.  Approved  for  public  release;  distribution  unlimited. 


Philosophical  Magazine 


3359 


displacement  field  of  the  screw  and  120°  dislocations,  as  compared  to  Figures  3 
and  4.  The  [111]  and  [112]  projections  are  shown  for  each  case.  In  Figure  5a,  the 
screw  dislocation  was  initially  positioned  a  [112]  unit  along  the  positive  ‘j’-axis 
(+0.8  nm),  away  from  the  origin.  In  Figure  5b,  the  screw  dislocation  was  initially 
positioned  a  [112]  unit  along  the  negative  ‘j’-axis  (— 0.8nm),  away  from  the  origin. 
In  both  cases,  after  relaxation,  the  screw  dislocation  is  partially  on  the  glide  (111) 
plane  and  partially  on  the  cross-slip  (111)  plane.  This  is  in  essence  the  presence  of 
cross-slip  nuclei  at  these  intersections.  These  results  are  similar  to  previous 
observations  of  athermal  cross-slip  nucleation  at  screw  dislocation  interactions 
with  Frank  interstitial  loops  in  fee  Ni  [18].  In  Figure  5a,  as  in  Figure  3,  the  screw 
dislocation  is  constricted  at  both  the  intersections  (negative  constrictions  [6,7]) 
around  one  of  the  Shockley  partials  of  the  120°  '/2[i01]  dislocation  residing  on 
the  (111)  plane.  Since  the  core  is  in  the  partially  cross-slipped  state,  a  companion 
positive  constriction  [6,7]  forms  between  the  intersections.  In  Figure  5b,  as  in 
Figure  4a,  the  a/6[Il2]  Shockley  partial  of  the  screw  dislocation  combines  with  the 
a/6[121]  Shockley  partial  of  the  120°  dislocation  to  form  an  extended  node  with  a 
resulting  short  Burgers  vector  of  the  type  1/6[011],  Lomer-Cottrell  barrier.  At  the 
other  end  of  the  Lomer-Cottrell  barrier,  the  screw  dislocation  is  constricted  to  form 
one  half  of  a  negative  constriction  and  is  cross-slipped  onto  the  glide  (111)  plane. 
Since  the  core  is  in  the  partially  cross-slipped  state,  a  companion  full  positive 
constriction  [6,7]  forms  between  the  intersections.  As  in  Figure  4a,  since  the 
interaction  between  the  1/6[011]  Burgers  vector  and  the  other  Shockley  partial  of 
the  120°  dislocation  is  weak,  the  120°  dislocation  needs  to  constrict  significantly 
to  attain  a  lower  separation  distance  at  the  extended  node  reaction  region.  Due  to 
the  fixed  boundary  conditions  along  ‘x’,  the  screw  dislocation  is  constricted  at  both 
ends  near  the  boundary  in  both  Figures  5a  and  5b,  forming  what  is  equivalent  to 
a  full  positive  constriction  [6,7].  In  this  case,  this  behavior  can  be  interpreted  as 
the  natural  formation  of  a  positive  constriction  inbetween  intersections  and 
no  energy  correction  is  required  to  obtain  the  energy  of  core  structures  PCSl  and 
PCS2. 


3.1.  Energetics  of  the  atomistic  intersection  core  structures 

The  energy  of  the  cells  was  evaluated  as  a  function  of  size  in  order  to  assess  the  errors 
introduced  by  the  finite  nature  of  the  simulation  volume  and  boundary  conditions. 
This  is  accomplished  by  calculating  the  energy  differences  between  core  structures 
truncated  at  different  finite  distances  along  the  ‘j’  and  ‘z’  directions.  Figure  6  shows 
a  plot  of  the  relative  energy  of  three  different  core  structures  obtained  for  the  screw 
—  120°  intersection  as  a  function  of  the  distance  along  ‘j’  and  ‘z’.  The  energy  is 
calculated  for  cells  of  lengths  5.0,  7.5,  10.0  and  12.5  nm  along  ‘j’  and  ‘z’  and  the 
energy  difference  obtained  between  the  different  core  structures  is  plotted  as  a 
function  of  these  cell  sizes.  Curve  (a)  in  Figure  6  is  a  plot  of  the  energy  differences 
between  two  identical  core  structures  obtained  for  the  intersection  (Figure  3, 
fully  residing  on  the  glide  (111)  plane)  with  two  different  initial  elastic  centers  and 
shows  that  the  measurement  of  relative  energy  has  an  error  of  the  order  of  0.2  eV 
using  this  technique.  The  partially  cross-slipped  configuration  (PCSl)  is  stable 
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Figure  6.  Relative  energies  of  the  core  structures,  Gil,  LC2  and  PCSl,  plotted  as  a  function 
of  the  truncation  distance  along  'y'  and  ‘z’.  Two  instantiations  of  the  core  structure  Gll  are 
also  compared  with  each  other  to  determine  the  accuracy  in  the  relative  energy  measurement. 
The  axes  dimensions  are  in  units  of  A. 


Table  3.  Relative  energies  of  various  core  structures 
obtained  atomistically  for  the  screw  dislocation  120° 
dislocation  intersection. 


Core 

Relative  energy  (eV) 

GLl 

0 

LCl 

5.3 

LC2 

0.8 

PCSl 

4.6 

PCS2 

5.8 

relative  to  the  cores  Gll  (fully  residing  on  the  glide  plane)  and  LC2  (fully  residing  on 
the  cross-slip  plane)  by  4-4.6  eV  and  the  energy  difference  values  are  fairly  converged 
at  a  distance  of  125  A  along  V  and  ‘z’.  If  one  subtracts  the  bulk  Stroh  constriction 
energy  (~1.9eV,  estimated  as  a  mean  of  the  negative  and  positive  constriction 
energies  at  a  stacking  fault  energy  of  90mJm~^  [6])  from  the  energy  of  the  core 
structures  Gll  and  LC2,  for  reasons  stated  in  Section  3.0,  one  obtains  that  the 
partially  cross-slipped  configuration  is  stable  relative  to  Gll  and  LC2  by  2. 1-2.7 eV, 
or  1-1.3  eV  per  intersection.  These  results  show  that  the  energy  of  a  negative 
constriction  is  significantly  lowered  at  the  intersections,  relative  to  the  bulk,  due 
to  the  stress  field  of  the  intersecting  dislocations.  Table  3  gives  the  relative  energy 
values  for  the  five  different  core  structures  shown  in  Figures  3-5  obtained  using 
the  technique  described  above  with  a  truncation  distance  along  ‘j’’  and  ‘z’  of  12.5  nm, 
without  the  Stroh  constriction  correction. 
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These  results  suggest  that  cross-slip  nuclei  are  readily  available  at  these 
intersections  and  that  cross-slip  nucleation  at  these  intersections  is  relatively  easy 
in  fee  Ni.  Also,  the  energy  of  core  structures  with  extended  nodes  (obtained  when  the 
screw  dislocation  was  initially  placed  a  [1 12]  unit  along  the  negative  ‘j’-axis,  LCl  and 
PCS2)  is  significantly  larger  than  the  core  structures  which  have  no  extended  nodes 
(obtained  when  the  screw  dislocation  was  initially  placed  a  [112]  unit  along  the 
positive  ‘j’-axis,  Gil,  LC2  and  PCSl).  This  may  be  related  to  the  constriction  of  the 
120°  dislocations  that  is  required  to  form  the  extended  node. 


3.2.  Check  of  initial  conditions  to  the  atomistic  simulations 

For  the  atomistic  simulations  described  above,  the  initial  separation  between  the 
screw  and  the  intersecting  120°  dislocations  was  kept  at  0.85  nm,  similar  to  what  was 
used  in  previously  published  atomistic  simulations  of  dislocation  intersections  [19]. 
In  order  to  assess  the  effect  of  the  initial  separation  distance  between  the  intersecting 
dislocations,  several  simulations  were  performed  with  increasing  initial  separation 
distance  between  the  intersecting  dislocations,  in  steps  of  0.85  nm,  up  to  a  separation 
distance  of  4.3  nm.  In  all  cases,  the  partially  cross-slipped  configuration,  with  the 
formation  of  cross-slip  nuclei,  was  observed,  even  when  the  stand-off  distance  exceed 
the  equilibrium  spacing  for  the  Shockley  partial  dislocations  of  a  relaxed  core. 
Figure  7  shows  a  plot  of  one  of  the  atomistic  core  structures  observed  for  the 
intersection  when  the  initial  separation  distance  between  the  intersecting  dislocations 
was  kept  at  3.4  nm.  This  suggests  that  a  small  initial  separation  distance  (~1.0nm) 
between  the  intersecting  dislocations  does  not  significantly  affect  the  atomistic 
simulation  results.  Investigation  of  further  increases  in  the  initial  separation  distance 
requires  larger  simulation  cell  sizes  than  used  here  because  of  the  fixed  boundary 
conditions  along  all  three  directions. 

Secondly,  the  initial  character  of  the  screw-character  dislocation  was  made  to 
deviate  from  the  perfect  screw  orientation  in  steps  of  2°  up  to  a  deviation  of  6°. 
Again,  in  all  cases,  the  partially  cross-slipped  configuration,  with  the  formation  of 
cross-slip  nuclei,  was  observed.  Figure  8  shows  a  plot  of  the  observation  of  a 
partially  cross-slipped  configuration  for  the  intersection  in  atomistic  simulations, 
when  the  initial  nearly  screw-character  dislocation  was  oriented  4°  away  from  the 
perfect  screw  orientation,  with  a  direction  of  [52  —48  —4].  Again,  further  tilting  of 
the  screw-character  dislocation  requires  larger  simulation  cell  sizes  than  used  here. 

Finally,  in  all  the  simulations  described  above,  two  120°  intersecting  dislocations 
having  a  separation  distance  of  30  nm  were  used,  in  order  to  minimize  the  effect 
of  fixed  boundaries  along  the  x-axis.  The  separation  distance  of  30  nm  corresponds 
approximately  to  a  dislocation  density  of  10*^  m“^,  which  is  representative  of  the 
upper  limits  of  densities  in  heavily  worked  materials  and  is  typical  of  what  may  be 
found  in  substructure  walls.  A  simulation  was  run  where  the  separation  distance 
between  the  two  intersecting  120°  dislocations  was  increased  by  a  factor  of  3, 
to  92.5  nm,  and  the  simulation  cell  size  along  the  .x-axis  was  also  increased  by  a  factor 
of  3,  to  185nm.  This  separation  distance  corresponds  approximately  to  a  dislocation 
density  of  lO'^^m”^.  Figure  9  shows  the  partially  cross-slipped  configuration  with 
the  formation  of  cross-slip  nuclei  observed  in  atomistic  simulations  under  these 
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Figure  7.  [1 1 1]  (x  plane)  and  [1 12]  (x  z  plane)  projections  of  the  partially  cross  slipped  core 
structure  for  a  screw  120°  intersection  in  fee  Ni  from  atomistic  simulations  when  the  initial 
separation  distance  between  the  screw  and  intersecting  dislocations  was  kept  at  3.4  nm. 
Atoms  with  eneigy  greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions 
are  in  units  of  A. 


conditions.  In  addition,  in  order  to  check  whether  the  presence  of  two  intersecting 
120°  dislocations  has  a  significant  effect  on  the  simulations  described  above,  a 
simulation  was  performed  with  only  a  single  intersecting  dislocation,  at  the  center  of 
the  simulation  cell.  Once  again,  the  partially  cross-slipped  configuration,  with  the 
formation  of  cross-slip  nuclei,  was  observed,  as  shown  in  Figure  10. 


3.3.  Periodic  boundary  condition  atomistic  simulations 

Several  simulations  were  performed  imposing  periodic  boundary  conditions  along 
the  screw  (x  or  110)  direction.  Such  boundary  conditions  could  be  imposed  only 
by  making  the  two  intersecting  dislocations  opposite  or  dipole  in  character  (-1-120° 
and  —60°).  Fixed  boundary  conditions  were  still  imposed  along  the  ‘j’  and  ‘z’ 
directions.  Figure  1 1  shows  a  plot  of  two  different  core  structures  obtained  for  the 
intersection  using  four  different  elastic  centers  for  the  initial  anisotropic  displace¬ 
ment  field,  with  periodic  boundary  conditions  imposed  along  the  ‘x’  direction. 
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4°  Inclination 
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Figure  8.  [11 1]  (x  plane)  and  [1 12]  (x  z  plane)  projections  of  the  partially  cross  slipped  core 
structure  for  a  screw  120°  intersection  in  fee  Ni  from  atomistic  simulations  when  the  initial 
direction  of  the  nearly  screw  character  dislocation  was  [52  48  4).  Atoms  with  energy  greater 
than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in  units  of  A. 


In  both  cases,  the  screw  dislocation  is  partially  cross-slipped  onto  the  cross-slip  (111) 
plane.  In  Figure  11a,  the  screw  dislocation  is  partially  cross-slipped  at  the  -f-120° 
intersection,  whereas  in  Figure  1  lb,  it  is  partially  cross-slipped  at  both  the  -j-]20°  and 
—60°  intersections.  These  results  suggest  that  the  partially  cross-slipped  configura¬ 
tion  forms  even  at  the  repulsive  —60°  intersection. 


3.4.  Dislocation  dynamics  simulations 

Figure  12  shows  a  plot  of  the  structure  of  a  120°  intersection  core  obtained  using  3D 
dislocation  dynamics  simulations.  The  dislocation  dynamics  simulation  results  are 
compared  with  the  corresponding  core  structure  obtained  for  the  120°  intersection 
with  atomistic  simulations.  Here,  the  screw  dislocation  was  initially  placed  in  the  3D 
dislocation  dynamics  simulation  cell  as  partially  cross-slipped  onto  the  cross-slip 
(111)  plane.  The  node  transitioning  the  screw  dislocation  from  the  (1 1 1)  glide  plane 
to  the  cross-slip  (111)  plane  was  free  to  move  along  the  screw  [110]  direction. 
Figure  12  shows  that  the  partially  cross-slipped  configuration  for  the  120° 
intersection  core  is  stable  and  good  correspondence  is  found  between  the  dislocation 
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Figure  9.  [1 1 1]  (jc  plane)  and  [1 12]  (jc  z  plane)  projections  of  the  partially  cross  slipped  core 
structure  for  a  screw  120°  intersection  in  fee  Ni  from  atomistic  simulations  when  only 
a  single  intersecting  dislocation  at  the  center  of  the  simulation  cell  was  used.  Atoms 
with  energy  greater  than  their  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions 
are  in  units  of  A. 


dynamics  and  atomistic  simulation  results.  We  conclude  that  the  elastic  interactions 
of  the  partial  dislocations  and  the  ensuing  recombination  of  the  partials  at  the 
junctions  capture  the  essential  details  of  the  dislocation  evolution  within 
the  atomistic  simulations.  Similar  dislocation  dynamics  simulations  were  found  to 
reproduce  all  the  five  core  structures  (Figures  3-5)  predicted  by  the  atomistic 
simulations. 


4.  Discussion 

The  results  of  this  work  clearly  show  that  several  dislocation  configurations  exist 
within  fee  crystals  containing  forest  dislocations,  where  a  glide  dislocation  can 
nucleate  a  cross-slip  event.  The  atomistic  simulation  results  show  that  cross-slip 
nuclei  are  readily  available  at  near-screw  dislocation  intersections  with  forest 
dislocations  residing  on  the  two  {111}  planes  other  than  the  cross-slip  plane  forming 
one  of  Lomer-Cottrell,  glide  or  Hirth  locks. 

This  mechanism  of  cross-slip  nucleation  is  in  addition  to  other  mechanisms 
of  obstacle  induced  cross-slip  initiation,  like  cross-shp  at  screw  dipoles,  jogs, 
surfaces,  etc.  [8,20,21].  However,  the  probability  of  dislocation  intersections  is 
relatively  large  as  compared  to  other  mechanisms  of  obstacle  induced  cross-slip. 
Also,  cross-shp  at  screw  dipoles,  jogs  and  surfaces  is  athermal  only  under  specific 
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[1-10] 

Figure  10.  [Ill]  (a:  y  plane)  and  [112]  (jc  z  plane)  projections  of  the  partially  cross  slipped 
core  structure  for  a  screw  120°  intersection  in  fee  Ni  from  atomistie  simulations  where 
the  separation  between  the  two  interseeting  120°  disloeations  was  92.5  nm.  Atoms  with 
energy  greater  than  the  energy  at  a  staeking  fault  are  plotted.  The  axes  dimensions  ate  in 
units  of  A. 


Figure  11.  [1 1 1]  (jc  y  plane)  and  [1 12]  (x  z  plane)  projeetions  of  two  different  core  structures 
obtained  for  screw  120°  and  60°  intersections  in  fee  Ni  from  atomistic  simulations  with 
periodic  boundary  conditions  along  the  screw  or  ‘.r’  direction.  Atoms  with  energy^  greater 
than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in  units  of  A. 
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Figure  12.  A  comparison  of  the  dislocation  dynamics  results  for  the  partially  cross  slipped 
core  structure  of  a  screw  120°  intersection,  with  atomistic  simulation  results.  [Ill]  (xy 
plane)  and  [112]  (x  z  plane)  projections  are  shown  for  each  case.  The  axes  dimensions  are 
in  units  of  Burgers  vector  magnitude  (2.5  A). 


conditions  [8,20,21].  This  is  an  important  result  since  it  qualitatively  justifies  the 
profuse  nature  of  cross-slip  in  fee  crystals,  unlike  the  current  models  for  thermally- 
activated  cross-slip  that  require  high  stress  to  provide  a  self-consistent  explanation 
[22].  This  finding  should  allow  higher-level  mesoscale  models  of  dislocation  behavior 
to  better  represent  the  cross-slip  process  without  resort  to  ad  hoc  assumptions 
about  obstacles. 

The  new  cross-slip  nucleation  model  has  a  variety  of  implications  for  crystal 
plasticity  in  fee  materials.  For  example,  within  the  present  mechanism,  the  frequency 
of  cross-slip  should  scale  with  the  forest  dislocation  density.  The  growth  of  such 
nuclei  should  depend  upon  the  relative  magnitude  of  local  stresses  on  the  glide  plane 
and  the  cross-slip  plane  at  the  partially  cross-slipped  screw  dislocation  intersection 
region.  Such  behavior  of  the  partially  cross-slipped  core  under  different  modes 
of  applied  stress  could  be  studied  using  atomistic  simulations  as  well  as  dislocation 
dynamics  simulations.  Still  it  would  be  instructive  to  perform  atomistic  simulations 
of  these  screw  dislocation  intersection  core  structures  in  other  fee  materials  with 
different  stacking  fault  energies  (i.e.  Al  and  Cu)  to  verify  the  results  in  Ni  as  well  as 
determine  the  energy  well  for  the  partially  cross-slipped  intersection  core  configu¬ 
ration  as  a  function  of  stacking  fault  energy.  The  stacking  fault  energy  may  also  play 
a  role  in  the  case  of  extending  these  configurations  under  stress  to  fully  accomplish 
the  cross-slip  process.  The  intersection  mechanism  of  cross-slip  nucleation  should 
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also  be  implemented  in  3D  dislocation  dynamics  simulations  as  an  alternative  to 
Escaig’s  model  in  dislocation-structure  evolution  methods  in  fee  materials.  Finally, 
it  would  be  instructive  to  revisit  the  Bonneville  and  Escaig  experimental  results  on 
cross-slip  in  fee  Cu  in  light  of  the  new  intersection  mechanism  for  cross-slip 
nucleation  [1,23]. 

According  to  Washburn  [10],  double  intersection  cross-slip,  where  the  segment 
that  has  been  pulled  into  the  cross-slip  plane  soon  encounters  another  attractive 
intersection  that  brings  it  back  onto  another  primary  glide  plane,  provides  a 
reasonable  mechanism  for  dislocation  multiplication  and  the  growth  of  slip  bands  at 
low  temperatures.  Naturally,  such  a  mechanistic  process  must  be  demonstrated 
via  modern  simulation  methods.  Also,  classical  theories  of  strain-hardening  assume 
that  dislocation  storage  in  Stage  II  of  single-crystal  fee  materials  is  a  result  of 
junction  formation  [24]  or  two-dimensional  (2D)  concave  loop  formation  [25]  as  the 
gliding  dislocation  traverses  through  an  array  of  forest  dislocation  obstacles  on  its 
glide  plane.  However,  one  of  the  major  problems  in  classical  strain-hardening 
models  is  to  explain  how  the  generation  of  a  3D  network  of  stored  dislocations 
occurs  as  a  consequence  of  2D  glide  [25].  We  note  that  the  intersection  cross-slip 
nucleation  mechanism  for  dislocation  storage  may  provide  a  convenient  mechanism 
of  generating  a  three-dimensional  network  of  stored  dislocations  from  2D  glide. 

The  anomalous  flow  behavior  of  intermetallics  (i.e.  Ni3Al)  has  thus  far  been 
modeled  using  core  transformations  that  require  thermal  activation.  Here  again, 
the  key  issue  has  been  that  constriction  energies  are  prohibitively  high  [3,4,26]. 
The  possibility  of  the  core  transformation  in  the  presence  of  interaction  with  other 
dislocations  provides  a  novel  mechanism  to  consider  and  explore  in  rationalizing  the 
flow  anomaly  of  these  types  of  materials.  However,  unlike  the  fee  metals,  it  should  be 
noted  that  to  rationalize  the  increase  in  flow  stress  with  temperature,  a  thermal 
component  to  the  locking  mechanism  is  required.  It  is  possible  that,  in  these 
materials,  the  dislocation  interaction  alone  is  insufficient  to  nucleate  core  trans¬ 
formation,  but  it  may  play  a  key  role  in  lowering  the  activation  energies  since  the 
explicit  formation  of  a  constriction  is  not  required. 

The  possibility  of  dislocation  interaction  induced  cross-slip  nucleation  raises 
the  question  of  why  work-hardening  is  found  to  saturate  with  plastic  strain.  The 
monotonic  increase  in  dislocation  density  with  plastic  strain  is  well  known  and  the 
attendant  increase  in  flow  stress  is  easily  rationalized.  However,  the  work-hardening 
rate  itself  decreases  with  strain.  If  dislocation  interactions  enhance  cross-slip 
nucleation,  how  can  a  decreasing  work-hardening  rate  with  increasing  dislocation 
density  be  rationalized?  One  possibility  is  that  there  is  a  critical  spacing  of  the  core¬ 
transforming  nuclei  beyond  which  there  is  no  enhancement  in  cross-slip.  This  spacing 
will  depend  on  how  fast  the  nucleus  grows,  which  in  turn  will  depend  on  the  stresses 
present  on  the  two  possible  glide  planes.  Both  atomistic  and  discrete  dislocation 
simulations  will  be  required  to  address  this  question  properly. 


5.  Summary 

These  atomistic  simulations  show  cross-slip  nucleation  can  be  easy  during 
dislocation-dislocation  interactions.  This  offers  a  pathway  for  using  atomistic 
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simulations  to  identify  types  of  interactions  that  result  in  cross-slip  and  use  these 
results  to  rationalize  in  higher  level  mesoscale  simulations  such  as  dislocation 
dynamics.  Summarizing  the  results:  for  a  120°  '/2[101]  Burgers  vector,  (Oil)  line 
direction  dislocation  residing  on  the  (111)  plane  intersecting  a  '/2[110]  screw 
dislocation,  the  following  observations  were  made  using  atomistic  simulations  of  the 
intersection  core  structure  using  embedded  atom  potentials: 

(1)  The  intersection  produces  five  different  possible  core  structures  where  the 
screw  dislocation  resides:  completely  on  the  glide  (111)  plane,  completely 
on  the  (111)  cross-slip  plane  in  two  core  structures  and  partially  on  both  the 
glide  (111)  and  cross-slip  (111)  planes  for  two  other  core  structures. 

(2)  One  of  the  partially  cross-slipped  configuration  for  the  core  structure  of  the 
intersection  is  stable  relative  to  the  other  two  core  structures  (residing  either 
fully  on  glide  (111)  or  fully  on  cross-slip  (111)  plane)  by  l.l-1.3eV  per 
intersection. 

(3)  The  extended  node  structure  for  the  intersections  is  relatively  unstable  with 
respect  to  the  constricted  structure,  presumably  due  to  the  formation  of 
significant  constriction  on  the  120°  intersecting  dislocation. 

(4)  Dislocation  dynamics  simulations  accounting  for  Shockley  partials  are 
shown  to  reproduce  the  products  of  dislocation  intersection  interacting  at 
the  atomistic  level. 
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